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I . SUMMARY 

^ The  objective  of  our  research  program  is  to  examine  the 

parameters  that  affect  the  seismic  signals  from  underground 
explosions.  Our  attention  is  primarily  directed  to  those 
features  of  the  seismic  waveforms  that  reliably  indicate  the 
explosion  yield.  Our  research  progrcun  includes  empirical 
studies  of  the  available  data,  experimental  studies  using 
small  charges  to  simulate  explosions  and  the  development  and 
application  of  theoretical  and  numerical  methods.  Emphasis 
is  on  the  latter,  in  particular  we  are  applying  techniques 
for  numerically  simulating  the  far-field  signals  from  both 
contained  and  cratering  underground  explosions.  The  numeri- 
cal simulation  techniques  represent  a synthesis  of  the  finite 
difference  methods  for  computing  ground  motion  in  the  near- 
source  large  displacement  regime  and  the  efficient  wave  prop- 
agation techniques  of  theoretical  seismology. 

During  the  second  three-month  period  of  our  present 
contract,  our  research  has  been  conducted  in  a number  of 
areas.  These  areas  and  the  major  results  obtained  in  each 
are  summarized  in  the  following  paragraphs. 

An  important  problem  engaging  our  attention  during  most 
of  this  quarter  was  an  attempt  to  understand  and  explain  the 
coupling  differences  between  explosions  in  difrerent  source 
materials  (granite,  salt,  saturated  tuffs  and  rhyolites) . 
Considerable  attention  was  given  to  the  effect  of  changing 
overburden  pressure  for  explosions  at  different  depths. 

Many  of  the  important  results  obtained  were  presented  orally 
at  several  DoD  meetiofc  during  the  quarter.  We  are  now 
preparing  a special  report  in  which  our  current  understanding 
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A number  of  research  efforts  were  brought  to  completion, 
or  nearly  so,  during  the  quarter.  This  report  is  devoted  to 
a summary  of  these  investigations.  We  reviewed  the  published 
evidence  bearing  on  the  nature  and  variadiility  of  the  anelas- 
tic  attenuating  properties  of  the  earth's  upper  mantle,  and  a 
topical  report  [Mitchell  and  Bache,  1977]  was  submitted.  We 
conclude  that  the  published  evidence  strongly  suggests  the 
presence  of  a low  velocity,  high  attenuation  zone  in  the  upper 
mantle  beneath  the  western  United  States.  This  zone  appears 
to  be  weak  or  absent  beneath  the  eastern  United  States  and 
central  Asia.  The  presence  of  such  a low  velocity  zone  is 
expected  to  significantly  lower  values  compared  to  values 
from  events  in  areas  where  it  is  absent.  However,  the  effect 
on  Mg  should  be  quite  small. 

A study  of  multiple  explosions  and  the  application 
of  certain  signal  analysis  techniques  to  separate  the  indi- 
vidual events  was  brought  to  completion  this  quarter  and  a 
topical  report  [Laimbert,  1977]  is  being  prepared.  The  object 
of  this  study  was  to  simulate  multiple  explosions  by  summing 
near-field  ground  motion  recordings  of  NTS  explosions  and 
then  attempt  to  identify  and  scale  the  individual  events 
by  analyzing  the  composite  signal.  Using  recordings  of  the 
large  Pahute  Mesa  events  MAST,  COLBY  and  POOL,  realistic 
multiple  explosion  signals  were  generated.  These  signals 
were  processed  by  our  program  MARS  (Multiple  Arrival  Recog- 
nition System)  which  uses  Hilbert  transform  envelopes  from 
narrow  band  filter  output  to  separate  individual  phase  arrivals 
in  a signal.  Using  close-in  instruments,  we  were  able  to 
identify  and  scale  the  individual  arrivals  down  to  quite 
small  time  separations.  Separation  times  can  be  resolved 
as  long  as  reasonable  signal  power  (signal/noise  > 1)  is  pre- 
sent at  frequencies  greater  than  3.5  divided  by  the  separation 
time.  For  exeunple,  if  good  signal  is  present  at  35  Hz,  phase 
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arrivals  separated  by  0.1  seconds  can  be  separated.  The 
results  studied  were  from  simulated  single  station  recordings 
of  a multiple  event.  For  an  actual  event,  it  would  be  desir- 
cOsle  to  employ  an  array  of  stations.  Analysis  of  such  data 
using  MARS  should  allow  identification  and  scaling  (by  yield) 
of  the  individual  events  in  the  multiple  explosion  array. 

We  have  been  conducting  an  experimental  program  in 
which  small-scale  modeling  of  underground  nuclear  explosions 
is  being  simulated.  The  tests  have  been  conducted  in  grout 
cylinders  and  the  explosives  are  0.25  grams  of  PETN,  which 
are  carefully  manufactured  to  be  spherically  symmetric  sources. 
In  our  last  quarterly  report  we  described  the  results  of  the 
first  set  of  experiments  in  which  the  explosions  were  fully 
contained.  A second  set  of  experiments  was  then  carried  out 
in  which  the  explosions  were  at  cratering  depths. 

One  important  reason  for  doing  the  experiments  described 
in  the  previous  paragraph  was  to  provide  a carefully  con- 
trolled data  set  for  testing  and  validating  our  numerical 
modeling  procedures.  We  have  modeled  the  contained  explo- 
sions using  a one-dimensional  finite  difference  code.  The 
code  and  the  constitutive  model  are  the  same  as  we  use  for 
modeling  underground  nuclear  explosions.  For  these  one- 
dimensional experiments  the  computed  "ground  motion"  time 
histories  agree  very  well  with  the  experimental  data.  This 
provides  important  validation  of  our  computational  procedures. 
We  are  now  engaged  in  modeling  the  cratering  experiments. 

Most  of  our  work  on  the  ground  motion  amplitude-yield 
scaling  for  underground  explosions  has  been  directed  toward 
explaining  the  teleseismic  (>2000  km)  data.  It  has  long  been 
known  that  the  increase  in  cunplitude  with  yield  is  markedly 
less  for  regional  (50-500  km)  body  wave  amplitude  than  for 
teleseismic  body  waves.  This  is  important  for  two  reasons. 
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First,  any  satisfactory  model  should  be  able  to  explain  all 
the  data.  Second,  the  regional  body  wave  data  sample  a 
different  portion  of  the  source  spectrum  than  other  data  and 
represent  another  opportunity  to  validate  calculations  of 
source  coupling.  For  these  reasons  we  have  made  a preliminary 
study  of  the  theoretical  amplitude-yield  scaling  for  the 
regional  data.  We  find  that  our  models  indeed  predict  a 
lower  slope  for  the  amplitude-yield  curve  in  this  region. 

The  primary  reason  is  that  the  regional  data  are  higher 
frequency  and  see  the  portion  of  the  source  spectrum  beyond 
the  corner  frequency.  This  effect  is  more  pronounced  at 
higher  yields  than  at  lower  yields,  thus  decreasing  the 
cimplitude-yield  slope. 
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II.  INTRODUCTION 

The  primary  objective  of  the  research  program  is  to 
systematically  examine  the  parameters  that  control  magnitude- 
yield  relationships  for  underground  nuclear  explosions. 

During  the  second  three-month  period  of  this  contract, 
research  has  been  conducted  in  the  following  areas: 

1.  A review  of  the  evidence  for  variability  in 
the  attenuating  properties  of  the  upper 
mantle  between  Eurasia  and  North  American 
(summarized  in  Section  3.1). 

2.  A study  of  the  utility  of  our  signal  analysis 
program  (MARS)  for  decomposing  and  scaling 
the  individual  events  in  a multiple  explo- 
sion array  (summarized  in  Section  3.2). 

3.  A laboratory  simulation  of  contained  and 
cratering  explosions  using  small  (0.25  gm) 
charges  in  grout  cylinders  (summarized  in 
Sec. ion  3.3) . 

4.  A finite  difference  numerical  modeling  of 
the  experiments  in  3 and  comparison  of 
observed  and  computed  data  (summarized  in 
Section  3.4). 

5.  A study  of  the  amplitude-yield  scaling  prop- 
erties for  large  underground  explosions 
recorded  at  regional  distances  (summarized 
in  Section  3.5). 

During  this  quarter  we  also  made  oral  presentations  at 
several  Department  of  Defense  (DoD)  meetings.  These  pre- 
sentations were  mostly  devoted  to  our  work  on  the  coupling 
of  explosions  in  various  source  materials  and  the  signifi- 
cance of  coupling  differences  for  mj^  and  M^.  This  work  is 
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central  to  our  contract  effort  and  most  of  our  attention  was 
devoted  to  these  issues.  A report  summarizing  the  current 
state  of  our  understanding  of  and  from  explosions  is 
now  in  preparation. 

The  remainder  of  the  report  is  devoted  to  technical 
discussions  of  results  in  the  five  areas  enumerated  above. 
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III.  TECHNICAL  DISCUSSION 

3.1  EVIDENCE  BEARING  ON  THE  NATURE  OF  THE  ANELASTIC  PROPER- 
TIES OF  THE  UPPER  MANTLE  BENEATH  NORTH  AMERICA  AND 

EURASIA 

During  this  quarter,  a topical  report  [Mitchell  and 
Bache,  1977]  was  prepared  sununarizing  work  relating  to  the 
seismic  attenuation  characteristics  of  the  western  United 
States  compared  with  those  of  the  eastern  U.S.  and  central 
Asia.  The  report  cites  previously  published  evidence  support- 
ing the  existence  of  a low  velocity,  high  attenuation  (low  Q) 
zone  in  the  upper-mantle  beneath  the  western  U.S.  (especially 
beneath  the  Basin  and  Range  Province)  and  the  lack  of  such  a 
zone  in  the  eastern  U.S.  and  central  Asia.  In  this  section 
we  summarize  the  main  conclusions  of  that  report. 

Models  of  seismic  velocity  and  Q inferred  from  data  for 
the  eastern  U.S.  require  a thin,  poorly  developed  low  velocity 
zone  or  none  at  all.  The  evidence  is  based  upon  long  seismic 
profiles  for  P waves  (reproduced  in  Figure  3.i  from  Masse 
[1973] ) , surface  wave  dispersion,  and  travel  time  residuals 
for  teleseismic  body  waves..  The  absence  of  a low  Q zone  is 
consistent  with  an  attenuation  study  done  for  the  upper  mantle 
beneath  the  eastern  U.S.  (see  Figure  3.2,  from  Lee  and  Solomon 
[1975] ) . 

Data  taken  in  the  western  U.S.  include  P and  S wave 
amplitudes  and  travel  time  residuals  from  teleseisms,  PL-P, 
travel  time  residuals,  regional  short-period  attenuation 
measurements  (from  the  nuclear  explosion  SALMON) , and  surface 
wave  attenuation  and  dispersion.  It  is  pointed  out  in  the 
report  that  in  North  America,  low  velocity  and  low  Q zones 
in  the  upper  mantle  coincide.  Also,  seismic  velocities  found 
in  the  low  velocity  zones  are  too  high  to  be  explained  by 
mineralogical  or  temperature  gradients.  The  low  velocities 
are  more  likely  due  to  partial  melting  or  water  in  the  upper 
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mantle,  both  of  which  result  in  higher  attenuation.  All  body 
wave  data  cited  (see  Figure  3.3)  support  the  existence  of  a 
low  velocity,  low  Q zone,  as  do  surface  wave  data  for  periods 
less  than  15  seconds.  Surface  wave  amplitudes  for  periods 
greater  than  15  seconds  do  not  exhibit  any  regional  variation, 
and  thus  measurements  should  not  vary  from  east  to  west  in 
the  U.S.  From  this  evidence  it  is  difficult  to  estimate  the 
extent  to  which  values  of  m^  will  vary,  but  low  Q zones  in- 
ferred are  sufficient  to  appreciably  lower  m^  values  in  the 
western  U.S.  It  is  also  expected  that  waves  entering  the 
western  U.S.  from  teleseismic  distances  would  be  further  at- 
tenuated elastically  by  the  seismic  velocity  contrast  between 
the  low  velocity  zone  and  its  upper  boundary.  Measurements  of 
high  heat  flow  and  electrical  conductivity  are  cited  as  further 
evidence  of  the  anomalous  nature  of  the  upper  mantle  beneath 
the  western  U.S. 

Although  data  for  central  Asia  and  northern  Europe  are 
sparse,  they  suggest  that  this  region  is  characteristic  of  a 
tectonically  stable  area.  Love  wave  group  velocities  were 
found  to  be  much  higher  than  in  the  neighboring  active 
Himalayas  and  Tibetan  Plateau.  Rayleigh  wave  attenuation  is 
similar  to  that  in  eastern  North  America.  These  and  other  data 
suggest  the  absence  of  a low  velocity,  low  Q zone  beneath 
Eurasia. 

It  is  concluded  that  the  evidence  indicates  the  pre- 
sence of  a low  velocity,  low  Q zone  in  the  western  U.S.  but 
not  in  the  eastern  U.S.  and  stable  areas  of  Eurasia  This 
zone  is  expected  to  affect  measurements  of  very  little,  but 
those  of  mj^  by  an  appreciable  amount. 

3.2  SIMULATION  AND  DECOMPOSITION  OF  MULTIPLE  EXPLOSION 

SEISMOGRAMS 

Lambert,  et  al . [1976]  and  Savino,  ^ al . [1976]  des- 
cribe the  results  of  some  experiments  on  the  separation  and 
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scaling  of  individual  events  from  seismic  records  of  multiple 
explosions.  Herein  we  summarize  the  continuation  of  that 
work.  A detailed  report  on  this  research  will  be  given  in  a 
topical  report  in  preparation.  The  purpose  of  the  study  was 
to  simulate  signals  from  multiple  explosions  and  to  develop 
techniques  to  identify  the  separate  events. 

The  techniques  used  to  decompose  signals  from  multiple 
events  are  based  on  applying  a sequence  of  narrow-band  filters 
over  a range  of  center  frequencies  to  these  signals  (see 
Lambert,  ^ al . , 1976,  for  details).  For  each  center  frequency, 
the  filtered  output  provides  a display  of  the  phase  (relative 
arrival  times)  and  amplitude  information  at  that  frequency. 
Resolution  of  smaller  time  difference  requires  high  center 
frequencies.  The  following  discussion  briefly  describes  two 
examples  of  our  most  recent  applications  of  these  techniques 
to  multiple  explosions. 


In  our  most  recent  work  we  simulated  signals  from  a 
multiple  event  for  which  all  signals  traversed  a common 
propagation  path  (Case  I)  and  for  which  the  propagation  paths 
were  different  (Case  II)  . Since  we  have  no  seismograms  for 
an  actual  multiple  event,  we  modeled  the  situation  in  Case  I 
by  using  records  of  a single  shot  at  a set  of  stations  as 
though  they  were  a set  of  shots  recorded  at  a single  station. 
That  is,  the  records  from  the  stations  were  simply  added 
together  to  form  a composite  seismogram,  which  was  then 
analyzed  to  extract  the  constituent  arrivals. 


The  events  used  for  Case  I were  the  Pahute  Mesa  events 
MAST  and  COLBY.  The  shot-station  geometry  is  depicted  in 
Figure  3.4.  Figure  3.5  shows  the  composite  signals  anc 
corresponding  Hilbert  transform  envelope  sums  (see  Lambert, 
et  al . , 1976)  for  various  combinations  of  records  for  the 
MAST  event.  The  arrival  times  at  the  respective  stations 
are  indicated  on  the  envelope  sum  plots.  The  signals 
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• Slapdown  / very  large  amplitudes 

• Largest  signal  amplitude  preceeding  slapdown 

O 2nd  largest  signal  amplitude  preceeding  slapdown 
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Figure  3.4, 


The  narrow-band  filter  output  for  the  multiple 
event  composed  from  MAST  records  SZ,  S2,  S3  and 
S4  is  displayed.  The  peaks  of  the  envelope  are 
plotted  versus  arrival  time  for  each  filter.  T 
eunplitudes  of  these  peaks  are  indicated  by  the 
symbols  described  in  the  legend. 


Aivplltude  (ft/sQC> 


Figure  3.5.  Composite  signals  (left  column)  and  Hilbert  trans- 
form envelopes  for  different  combinations  of  re- 
cording stations.  The  onset  time  at  the  stations 
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corresponding  to  the  arrivals  at  SZ/S2,  S3  and  S4  are  seen  as 
distinct  bumps  on  at  least  one  of  the  envelope  sums.  The 
difference  in  arrival  times  between  SZ  and  S2  is  too  small 
(0.05  second)  to  be  resolved  on  this  plot.  The  large  spikes 
at  the  S5  arrival  time  are  due  to  spall  closure  and  obscure 
the  S5  arrival.  To  further  decompose  the  composite  seismo- 
grams, we  refer  to  the  narrow-band  filter  outputs  shown  in 
Figure  3.6.  A coherent  event  is  identified  as  a peak 
(identified  by  various  symbols  on  the  figure)  which  persists 
at  several  adjacent  frequencies  and  times.  For  example,  the 
event  labeled  "largest  amplitude  preceding  slap-down"  shows 
at  the  same  time  for  each  frequency.  Note  that  for  high 
frequencies,  an  event  preceding  the  largest  before  slap-down 
emerges.  This  has  the  travel  time  for  SZ  and  we  see  that 
the  technique  has  resolved  the  S2  and  SZ  arrivals.  In  general, 
our  results  show  that  separation  of  separate  arrivals  becomes 
possible  when  the  narrow-band  center  frequencies  are  above 
about  3.5/T,  where  T is  the  difference  in  arrival  times. 

Various  other  coherent  arrivals  may  be  seen  on  this  plot  that 
may  be  attributed  to  propagation  path  effects.  Other  inter- 
pretations and  analyses  for  different  composite  seismograms 
are  discussed  in  the  topical  report  now  in  preparation. 

We  proceed  now  to  Case  II,  where  seismograms  from  several 
events  were  used  to  make  up  simulated  multiple  explosion  records. 
We  used  the  events  MAST,  COLBY  and  POOL.  The  stations  used 
for  COLBY  and  POOL  were  on  a line  between  the  shot  points  and 
thus  are  likely  to  have  aspects  of  their  travel  path  in  common. 
The  travel  path  to  MAST,  on  the  other  hand,  is  likely  to  be 
rather  different.  Figure  3.7  shows  the  geometry  for  this 
problem.  The  stations  are  positioned  in  a line  according 
to  the  distance  from  their  respective  sources.  As  in  Case  I, 
we  plot  the  composite  and  envelope  sums  (Figure  3.8)  with 
the  arrival  times  indicated  on  the  plot.  One  can  see  distinct 
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Figure  3.7.  Vertical  section  of  geometry  for  Case  II. 
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peaks  (not  evident  on  the  composite  seismograms)  on  the 
envelope  sum  plot  for  COLBY  S2  and  S4  and  MAST  S3,  with  a 
much  less  distinctive  bump  for  POOL  S9.  Thus,  for  the  second 
case,  the  technique  proved  successful  in  separating  the 
arrivals . 

These  examples  and  several  others  are  described  in  our 
forthcoming  topical  report.  We  find  that  the  narrow-band 
filtering  technique  would  be  of  great  value  for  identifying 
and  scaling  the  individual  events  in  a salvo  of  explosions. 

3 . 3 A LABORATORY  SIMULATION  OF  CRATERING  EXPLOSIONS 

In  our  previous  quarterly  technical  report  [Bache, 
et  al . , 1977] , we  described  some  small-scale  modeling  experi- 
ments designed  to  provide  benchmark  data  for  computer  simu- 
lations of  underground  explosions.  In  those  experiments,  we 
obtained  displacement-time  histories  on  a free  surface  due  to 
a small  spherical  charge  detonated  within  a uniform  medium  of 
known  properties.  Our  first  results  were  for  fully-contained 
explosions.  In  our  most  recent  tests,  the  charge's  depth  of 
burial  was  reduced  and  cratering  occurred. 

Figure  3,9  shows  the  test  geometry.  The  medium,  a 
high-strength  concrete,  is  poured  in  the  form  of  a right 
circular  cylinder  of  diameter  122  cm.  Over  one  end  of  the 
cylinder,  the  motion  transducers  were  arrayed  at  the  vertices 
of  a 30  cm  equilateral  triangle.  About  31  cm  below  each  vertex 
of  the  triangle,  a spherical  charge  (1/4  gram  of  PETN)  was 
buried.  The  charges  were  fired  one  at  a time,  thus  each  con- 
crete cylinder-transducer  setup  yielded  three  complementary 
sets  of  data. 

Table  3.1  summarizes  the  experiment  matrix.  Figures 
3,5.1  and  3.5.2  of  our  previous  report  show  typical  results 


Vertical 
direction  is 
Z axis. 


Vertex  1 
Vertex  2 
Vertex  3 

VI,  Vertical  Displacement  Gauge 
HI,  Horizontal  Displacement  Gauge 

V2 , Vertical  Displacement  Gauge 
A2,  Vertical  Accelerometer 

A3 , Vertical  Accelerometer 

Charge  1 
Charge  2 
Charge  3 

Cylinder  Bottom 


(-15.2,  8.7,  0) 

(15.2,  8.7,  0) 

(0,  -17.6,  0) 

(-14.0,  8.7,  0) 

(-16.5,  8.7,  -0.5) 
Sensitive  axis  I lx 
(14.0,  8.7,  0) 

(16.5,  8.7,  0) 

(0,  -17.6,  0) 

(-15.2,  8.7,  -31) 
(15.2,  8.7,  -31) 

(0,  -17.6,  -31) 

(-,  -,  -33  or  -37  or 
-60  [Table  3.1]) 


Figure  3.9.  Concrete  test  bed. 
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for  the  contained  shots  and  Figures  3.10-3.12*  of  the  present 
report  give  the  results  for  the  cratering  events.  In  general, 
there  is  good  agreement  in  the  wave  shapes  between  the  three 
shots  in  each  configuration,  confined  or  cratering.  However, 
there  are  differences  in  absolute  amplitudes  which  are  larger 
than  our  estimated  errors  of  ±0.3  ym.  For  the  confined  shots, 
we  believe  these  differences  are  due  to  slight  variations  in 
the  coupling  between  the  concrete  and  the  charge  surface. 

For  the  subsequent  cratering  tests,  our  emplacement  technique 
for  the  charges  was  improved,  and  it  does  appear  that  the 
peak  displacements  are  correlated  with  the  crater  size  and 
depth  of  burial.  Figure  3.14  shows  the  measured  displacement 
time  history  from  a contained  shot.  This  displacement  gauge 
was  located  30.5  cm  from  the  explosive  source.  The  gauges 
labeled  VI  (Figure  3.10)  and  V2  (Figure  3.11)  measured  free 
field  displacement  ground  motion  from  an  explosive  source 
which  produced  a crater.  The  gauges  were  located  31.3  cm  and 
31.6  cm  from  the  explosive. 

The  major  difference  in  the  radiated  ground  motion  be- 
tween the  contained  and  cratering  modes  (Figures  3.10,  3.11  and 
3.14)  is  the  large  rebound  produced  by  the  cratering  shots. 

This  is  probably  caused  by  the  release  of  the  highly  stressed 
region  around  the  cavity  by  the  free  surface  reflection.  It 
is  very  important  that  this  release  mechanism  be  understood  in 
terms  of  its  effect  on  and  M^. 

3.4  NUMERICAL  SIMULATION  OF  MODEL  EXPERIMENTS 

In  our  preceding  quarterly  report  [Bache,  e;t  al . , 1977], 
we  described  small-scale  experiments  which  modeled  underground 

*The  error  bars  shown  include  the  estimated  uncertainties  due 
to  gauge  and  recorder  calibrations  and  electrical  noise.  The 
figures  also  quote  the  times  of  first  arrival  and  zero  crossing 
and  first  peaks  in  the  accelerometer  records.  Integrations  of 
A2 ' s data  are  not  available  because  of  limitations  in  our 
recording  apparatus. 
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contained  explosions  and  included  measurements  of  displace- 
ment-time histories  at  the  free  surface.  These  experiments 
have  now  been  numerically  simulated  using  material  properties 
from  laboratory  measurements  and  material  models  most  of 
which  were  described  in  an  earlier  report  [Bache,  e^t  al  ■ , 1975] 
One  change  in  the  source  calculation  method  was  a new  effec- 
tive stress  model  which  does  not  depend  on  a relaxation  time. 
The  elimination  of  this  time  parameter  is  important  for  the 
scaling  of  source  calculations. 

The  free  surface  ground  motion  predicted  by  the  numeri- 
cal simulation  of  the  model  experiments  (shown  in  Figure  3.13) 
gives  excellent  agreement  with  the  measured  displacement-time 
histories  (Figure  3.14).  Only  one  calculation  was  done  so 
this  agreement  was  not  achieved  by  a posteriori  adjustment  of 
material  models  or  material  properties.  It  therefore  pro- 
vides important  confirmation  of  our  source  calculation 
methods.  A more  complete  description  of  this  numerical  simu- 
lation is  given  in  the  following  paragraphs. 

CONFIGURATION 

The  model  experiments  consisted  of  a spherical  1/4-gram 
charge  of  powdered  PETN  loosely-pacJced  to  a density  of  1 gm/cm^ 
The  containing  material  was  "high  early"  grout  formed  into  a 
right  circular  cylinder  with  a diameter  of  120  cm  and  thick- 
ness of  60  cm.  The  motion  was  measured  at  one  end  of  the 
cylinder  and  the  charges  were  located  30.5  cm  below  the  sur- 
face. Numerically  this  configuration  was  modeled  as  shown 
in  Figure  3.15  with  a spherical  charge  in  a sphere  of  concrete. 

GROUT  MODELING 

The  major  types  of  material  behavior  which  we  have 
included  in  modeling  the  grout  are  those  attributable  to 
water  content,  air-filled  porosity,  brittle  failure  and  plastic 
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flow.  The  methods  for  simulating  these  phenomena  are  the 
same  as  those  described  previously  [Bache,  e^  al . , 1975] 
except  for  the  treatment  of  effective  stress  described  below. 

The  basis  of  an  effective  stress  law  is  the  assumption 
that  the  fluid  pressure  in  fluid  saturated  rocks  should  be 
included  in  the  calculation  of  the  stress  state  which  deter- 
mines material  strength.  The  definition  of  stress  state,  P, 
which  we  have  adopted  is 

^ p - H-t)  (3.1 

where  P is  the  pressure  component  of  the  stress  tensor  and  J' 
is  the  third  deviatoric  stress  invariant 


J3  = {0^^  + P)  (0^2  + P)  (033  + P)  . (3.2 

Here,  °22'  *^33  principal  stresses.  For  our  new 

effective  stress  method  the  stress  state  P is  adjusted  by 
subtracting  the  pore  fluid  pressure,  P^,  before  the  yield 
function  is  evaluated.  The  material  strength  is  determined 
by  the  following  yield  function: 


Y P ■ Pf 

Y = Yf.  + -^  (P  - P„)  2 

Up  r p 

m m 


where  P , Y and  Y„  are  material  constants, 
mm  0 

The  fluid  pressure  P^  is  evaluated  as  the  following 
function  of  the  crack  strains,  "^22'  ^33'  the  material 

distension,  a,  which  is  the  ratio  of  volume  of  the  solid 
material  to  volume  of  material,  including  air-filled  voids: 
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oi,p  — ot  ( 1 + + ^22  ^33^  ' 

= P for  = 1 

(O^m  “ a \ 

-jr^j  1 ^ “t  “e 

P-  = 0 for  a_  > a 
f T e 


(3.4) 


In  this  function,  P^  = P for  a,j,  = 1 because  the  fluid  pressure 
is  in  equilibrium  with  the  rock  constituent  when  all  air- 
filled  porosity  is  removed  from  the  material.  For  between 
1 and  the  elastic  limit,  a^,  some  of  the  air-filled  porosity 
has  been  irreversibly  removed  and  P^  is  taken  as  a linear 
function  of  a_.  For  a_  > a the  fluid  pressure  is  zero  because 
the  pores  are  elastically  deformed  but  not  irreversibly  crushed. 

The  parameters  defining  the  material  behavior  for  grout 
for  this  simulation  are  summarized  in  Table  3.2.  These  param- 
eters were  derived  from  our  measurements  of  the  P-wave  velocity 
Vp,  density  p and  the  unconfined  compressive  strength  Y^. 
Additional  characteristics  were  obtained  from  the  extensive 
data  available  on  DF-5a  grout,  [Cherry  and  Peterson,  1970] , 
which  is  a similar  material.  In  Table  3.1,  the  parameters  A 
and  B define  the  equation  of  state  in  the  form 


Au  + Bu 


where  y = — - 1.  The  values  for  shear  modulus,  G,  and  for 

*^0 

bulk  modulus.  A,  were  determined  from  V and  p by  assumina  a 

p o 

Poisson's  ratio  of  0.1,  which  is  appropriate  for  this  type  of 
material . 


EXPLOSIVE  SOURCE 

The  explosive  source  was  PETN  powder  loosely-packed  to 
a density  of  1 gm/cm^.  (Solid  PETN  has  an  initial  density  of 
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TABLE  3.2  MATERIAL  PARAMETERS  USED  FOR 
HIGH  EARLY  GROUT 


u 

A 

B 


m 


m 

G 


2.52  gm/cm^ 

0.345  kbar 
3.38  X 10^  cm/sec 
106.5  kbar 
67.3  kbar 
0.1725  kbar 
3 kbar 
5 kbar 
137  kbar 


P = 10  kbar 
c 

P =0.04  kbar 
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1.77  gm/cm^  and  a Chapman-Jouguet  density  of  2.40  gm/cm^.) 
The  results  of  Finger,  et  al . [1976] , provide  the  values 
for  the  constants  in  the  JWL  equation  of  state  for  the  pro- 
ducts of  detonation  of  PETN, 


P = A"  (1  - 


-R- V 

) e ^ (1 


R,V 


ti)  . _ 2 , toE 

r::v^  ^ ^ ^ ' 


(3.5) 


and  for  their  isentropic  expansion, 


P 

s 


A'e 


-R^V 


+ B'e 


-R2V 


V 


,(i)+l 


(3.6) 


In  Equations  (3.5)  and  (3.6),  V is  the  specific  volume,  1/p, 
and  the  constants  A^,  Q',  C',  R^^,  R2 , and  uj  are  given  in 
Table  3.3. 

For  the  numerical  simulation,  the  source  was  modeled  by 
a pressure  loading  function  defined  by  Equation  (3.6)  and 
applied  at  the  cavity  boundary.  The  pressure  function  was 
started  with  the  PETN  detonation  products  expanded  to  the 
actual  volume  of  the  loose-packed  PETN  which  corresponds  to  an 
initial  radius,  R^,  of  0.39  cm.  Thus,  the  initial  pressure 
was  31.2  kbar,  which  is  well  below  the  Chapman-Jouguet  pres- 
sure of  320  kbars. 


COMPARISON  OF  THE  RESULTS 

The  agreement  between  the  measured  and  calculational 
results  shown  in  Figures  3.13  and  3.14  is  considered  to  be 
very  good.  No  adjustment  of  models  or  material  properties 
was  necessary.  The  difference  between  the  calculated  dis- 
placement and  the  measured  displacement  is  within  the  range  of 
accuracy  of  the  measurements.  This  is  demonstrated  by  the 
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TABLE  3 . 3 JWL  CONSTANTS  FOR  PETN 


A' 

— 

7960 

kbar 

B' 

= 

192  : 

kbar 

C' 

= 

6.65 

kbar 

= 

8.50 

gm/cm 

= 

2.12 

gm/cm 

(ii 

= 

0.25 

^cj 

= 

320  : 

kbar 

^cj 

= 

2.40 

gm/cm 

= 

1.77 

gm/cm 
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data  from  a second  experiment  plotted  in  Figure  3.16.  The 
displacement  time  histories  shown  are  from  two  displacement 
gauges  and  a twice-integrated  accelerometer  record.  The 
three  gauges  were  in  close  proximity  and  should  be  measuring 
the  same  motion.  The  test  configuration  for  these  data  was 
the  same  as  that  for  the  data  of  Figure  3.14  so,  from  a 
theoretical  point  of  view,  the  displacements  in  Figures  3.14 
and  3.16  should  be  identical.  However,  the  peak  displacements 
range  between  4.5  um  and  6.7  utn  while  the  rise  time  varies 
between  14  and  20  tsec.  Our  calculated  displacement-time 
history  shows  a peak  displacement  of  6.6  um  and  a rise  time 
of  20  usee. 

3.5  REGIONAL  SCALING  OF  SEISMIC  AMPLITUDE  VERSUS  YIELD 

The  purpose  of  this  section  is  to  describe  some  pre- 
liminary results  for  research  on  scaling  of  seismic  amplitudes 
with  yield  for  nuclear  explosions  recorded  at  regional  dis- 
tances (100  to  1000  km) . It  has  been  pointed  out  by  Springer 
and  Hannon  [1973] , that  the  slope  of  log  amplitude-log  yield 
plots  is  not  unity  for  Pahute  Mesa  explosions  recorded  at 
regional  ranges.  Figure  3.17  shows  a plot  of  "b"  amplitude 
versus  yield  for  Pahute  Mesa  explosions  recorded  at  MINA, 
Nevada.  The  b amplitude  is  measured  from  the  first  cycle 
on  the  record  as  shown. 


The  values  on  the  plot  (provided  by  Lawrence  Livermore  Labora- 
tory) are  amplitudes  corrected  by  the  aain  of  the  instrument 
at  1 Hz,  and  have  not  been  corrected  for  the  response  of  the 
instrument  as  a function  of  frequency.  A straight  line  indicat- 
ing the  trend  shows  that  the  slopes  on  the  log-log  plots  are 
indeed  less  than  one  ("^  0.7). 
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Figure  3.16, 


Free  surface  vertical  displacement  measured  in 
three  different  ways.  The  largest  peak  is  from 
a twice- integrated  accelerogram,  while  the 
other  two  displacements  are  from  displacement 
gauges . 
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Towards  understanding  this  data,  the  following  cal- 
culations were  done.  First,  the  sources  for  Pahute  Mesa 
were  represented  by  reduced  displacement  potentials  (RDP) 
found  numerically  and  used  in  previous  teleseismic  amplitude 
studies  [see,  e.g.,  Bache,  et  al . , 1976].  The  RDP's  at 
different  yields  were  found  by  cube  root  scaling.  The  RDP’s 
(shown  in  Figure  3.18)  were  propagated  through  a simple 
earth  structure  (Figure  3.19)  to  a range  of  200  )cm  using  a 
generalized  ray  program.  The  results  are  then  convolved  with 
an  LRSM-SP  instrument  response.  Figure  3.20  shows  the  com- 
parisons of  the  theoretical  results  (using  source  204)  with 
observations.  The  slopes  are  seen  to  be  in  good  agreement. 

We  shall  discuss  these  comparisons  further,  but  shall  first 
discuss  the  factors  controlling  the  slope  of  the  log  b-log  Y 
curves  (hereinafter  referred  to  as  versus  Y curves) . 

Several  effects  contribute  to  the  slopes  of  the  b 
versus  Y curves.  The  first  is  the  filtering  of  the  signal 
by  the  instrument.  The  short-period  instrument  behaves  as  a 
narrow-band  filter  centered  about  3 Hz.  As  the  yield  changes, 
the  instrument  samples  different  parts  of  the  source  spectrum. 
This  could  be  visualized  by  thinking  of  the  instrument  response 
as  fixed,  and  sliding  the  source  spectrum  to  the  left  as  yield 
increases.  Thus,  the  second  feature  controlling  the  slope  of 
b versus  Y is  the  shape  of  the  source  spectrum.  To  see  this, 
consider  Figure  3.21  in  which  b versus  Y curves  are  plotted 
for  RDP's  of  different  shape,  one  being  strongly  peaked  near 
the  corner  frequency  and  one  fairly  flat.  The  effects  of  the 
crustal  structure  are  not  included  and  the  anelastic  attenuation 
factor  is  T/Q  = 0.1  for  both  plots.  The  curves  become  parallel 
at  low  or  high  yields  but  differ  markedly  for  intermediate 
yields  due  to  the  different  shape  of  the  source  functions 
near  the  corner  frequency. 
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Figure  3.18.  Amplitude  spectra  of  reduced  displacement  poten- 
tials for  three  materials,  the  amplitude  is 
scaled  to  0.02  kt  while  the  frequencies  have  been 
scaled  to  10  kt. 
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200  km 


Figure  3.19.  Schematic  of  ray  paths  and  source-receiver 
geometry.  Wavy  line  indicates  a head  wave. 
Vertical  traverses  are  P waves.  Values  in 
parentheses  are  the  P and  S wave  velocity 
and  density. 
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Another  factor  controlling  the  slope  of  b versys  Y 
curves  is  the  anelastic  attenuation.  Larger  values  of  T/Q 
result  in  greater  attenuation  of  high  frequencies  relative 
to  low  frequencies.  Thus,  as  T/Q  increases,  amplitudes  from 
lower  yield  events  are  attenuated  relative  to  higher  yield 
events,  which  results  in  steeper  b versus  Y curves.  Figure 
3.22  illustrates  the  effect  of  changing  T/Q  from  0.5  to  0.1 
for  source  254  with  the  crustal  response  removed. 

The  final  feature  we  shall  discuss  is  the  effect  of 
the  crustal  response.  For  our  simple  model  at  the  range 
200  km,  the  response  is  composed  of  the  head  wave  from  the 
direct  P from  the  source  and  the  head  wave  from  the  free  sur- 
face reflected  (pP)  phase  (see  Figure  3.19).  The  signal 
radiated  from  a head  wave  at  an  interface  has  a time  history 
which  is  essentially  the  integral  of  the  incident  wave.  This 
is  shown  in  Figure  3.23,  where  b versus  Y curves  are  plotted, 
for  source  254,  with  T/Q  = 0.1,  with  and  without  the  crustal 
response.  It  can  be  seen  that  the  slopes  are  reduced  by  the 
crustal  response. 

We  now  return  to  our  comparison  of  the  theoretical  and 
observed  curves  in  Figure  3.20.  As  described  above,  our 
theoretical  amplitude-yield  slope  is  in  close  agreement  with 
that  of  the  observations,  indicating  that  the  main  features 
controlling  the  slope  have  been  included  in  the  theoretical 
model.  The  theoretical  amplitudes  are,  however,  low  by  a 
factor  of  2 or  3.  Considering  that  no  attempt  was  made  to 
refine  crustal  models  and  values  of  T/Q  beyond  the  first-cut 
models  discussed  here,  nor  to  determine  the  station  character- 
istics of  MINA,  nor  to  include  effects  of  depth  of  burial 
on  source  material  properties  (only  cube  root  scaling  was 
used) , this  is  a reasonable  discrepancy.  We  point  out  that 
a number  of  factors  that  influence  the  amplitudes  have  not 
been  included  as  yet.  The  crustal  model  is  extremely  simple 
and  only  a "first-cut"  estimate  of  Q was  used.  We  will  do  some 


43 


(10  11) 


2 5 10  20  50  100  200  500 

Yield  (kt) 


Figure  3.23.  Log  b versus  log  yield  for  source  254,  T/Q  = 
0.1,  with  (top)  and  without  (bottom)  the  head 
wave  crustal  response.  The  curves  have  been 
shifted  vertically  for  ease  of  comparison. 
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